An effective chiral Hadron-Quark Equation of State 
Part I: Zero baryochemical potential 
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We construct an efTective model for the QCD equation of state, taking into account chiral sym- 
metry restoration as well as the deconfinement phase transition. The correct asymptotic degrees of 
freedom at the high and low temperature limits are included (quarks ^ hadrons). The model shows 
a rapid crossover for both order parameters, as is expected from lattice calculations. We then in- 
vestigate the thermodynamic properties of the model at = 0. All thermodynamic quantities are 
in qualitative agreement with lattice data, while apparent quantitative differences can be attributed 
to hadronic contributions and excluded volume corrections. 



INTRODUCTION 

The recent experimental results at the Relativistic 
Heavy Ion Colhder (RHIC), suggesting the creation of a 
nearly perfect fluid [li 0, S Q , have fueled interest in the 
study of bulk properties of strongly interacting matter 
(QCD). Heavy ion experiments at different beam ener- 
gies try to map out the QCD phase diagram, especially 
the region where one expects a phase transition from a 
confined gas of hadrons to a deconfined state of quarks 
and gluons (QGP) 0, S 0, S, i, Ei El, [3, 0, [3, [H . 
To relate any experimental observables to the properties 
of the matter produced in heavy ion collisions, a pro- 
found understanding of the thermodynamics of QCD has 
to be obtained and integrated in model simulations of 
these collisions. 

One can define two different phase transitions, the first 
being the chiral phase transition associated with chiral 
symmetry restoration in the vanishing quark mass limit, 
where the chiral condensate serves as a well defined order 
parameter. In the limit of heavy quarks, a deconfinement 
phase transition with the Polyakov loop as order parame- 
ter, is assumed. Physical quarks however have intermedi- 
ate masses and one would expect that at least the decon- 
finement order parameter is not so well defined anymore. 
There could even be some mixing of the two order param- 
eters, accounting for some lattice QCD observation that 
deconfinement and chiral restoration occur at the same 
temperature (at least at = 0) [13, H H 0, E | . 
Such lattice calculations at finite temperature are an im- 
portant tool for the investigation of the QCD phase dia- 
gram. For the thermodynamics of the pure gauge theory 
high accuracy data is available 22], and the equation of 
state (EoS) of strongly interacting matter at vanishing 
chemical potential is reasonably well understood 2^, 23| ■ 
Here lattice predicts a rapid crossover for the deconfining 
and chiral phase transitions. 

At finite baryo-chemical potential, lattice calculations 



suffer from the so called sign problem. There are sev- 
eral different approaches to obtain results at finite /is 



;rai ainerent approacnes to obtain results at nnite /is 
U, m, S [13,121, [H, [13, but yet no clear picture, 
especially about the existence and location of a possible 
critical end point, has emerged. 

Recent considerations based on connecting the large Nc 
limit with real-word QCD draw an even more exotic pic- 
ture of the phase diagram, where the critical tempera- 
tures of the deconfinement and chiral phase transitions 
disconnect and depart in the region of high net baryon 
densities ^32l] . 

In our approach we combine, in a single model, a well- 
established flavor-SU(3) hadronic model with a quark- 
gluon description of the highly excited matter. This al- 
lows us to study the chiral-symmetry and confinement- 
deconfinement phase structure of the strongly interacting 
matter at high temperatures and densities. In addition 
we obtain an equation of state of hadronic and quark 
matter that is applicable over a wide range of thermo- 
dynamical conditions and that can therefore be used in 
heavy- ion simulations with very different beam energies. 



MODEL DESCRIPTION 

In our approach we derive the EoS of hot and dense 
nuclear matter using a single model for the hadronic and 
quark phase. The model includes the correct asymptotic 
degrees of freedom, namely a free gas of quarks and 
gluons at infinite temperature, and a gas of hadrons 
having the correct vacuum properties at vanishing 
temperature. The model also predicts the structure 
of finite nuclei, nuclear and neutron matter properties 
and a first order liquid-vapor phase transition. The 
two phase transitions that are expected from QCD, the 
chiral and deconfinement transitions, are also included 
in a consistent manner. 

In the following we will show how we describe the 
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different phases of QCD and how we combine them in a 
single model. 

We describe the hadronic part of the EoS, using a 
flavor-SU(3) model which is an extension of a non-linear 
representation of a sigma-omega model including the 
pseudo-scalar and vector octets of mesons and the bary- 
onic octet and decuplet (for a detailed discussion see 

The Lagrangian density of the model in mean field ap- 
proximations reads: 



L 



(1) 



where besides the kinetic energy term for hadrons, the 
terms: 



Li 



■EiV'i[7o(5»(^w-}-5,;0(/) + m*]?/;j, (2) 



masses of the baryons, nuclear saturation properties and 
asymmetry energy as well as the A-hyperon potentials. 
The vacuum expectation values of the scalar mesons are 
constrained by reproducing the pion and kaon decay 
constants. 

The extension of the hadronic SU(3) model to quark 
degrees of freedom is constructed in analogy to the PNJL 
model [371, l38[ . The sigma model uses the Polyakov loop 



$ as the order parameter for deconfinement. $ is defined 
via $ = iTr[exp {i J dTA4)], where A4 = iAq is the tem- 
poral component of the SU(3) gauge field. One should 
note that one must distinguish and its conjugate <&* 
at finite baryon densities [ill, [sl, [3] , as they couple dif- 
ferently to the quarks, respective antiquarks. 
In our approach the effective masses of the quarks are 
generated by the scalar mesons except for a small ex- 
plicit mass term {Snig — 5 MeV and Siris — 105 MeV for 
the strange quark): 



Xa 



(3) 



represent the interactions between baryons and vector 
and scalar mesons, the self-interactions of scalar and vec- 
tor mesons, and an explicitly chiral symmetry breaking 
term. The index i denotes the baryon octet and decuplet. 
Here, the mesonic condensates (determined in mean-field 
approximation) included are the vector-isoscalars lu and 
(/), and the scalar-isoscalars a and ^ (strange quark- 
antiquark state). Assuming isospin symmetric matter, 
we can neglect the p-meson contribution in Eq. [3l 

The last four terms of JSl were introduced to model 
the QCD trace anomaly [34], where the dilaton field x 
can be identified with the giuon condensate. 

The effective masses of the baryons (of the octet) are 
generated by the scalar mesons except for an explicit 
mass term {Sm^ = 120 MeV): 



(4) 



while, for simplicity and in order to reduce the number 
of free parameters, the masses of the decuplet baryons 
are kept at their vacuum expectation values. With the 
increase of temperature/density, the a field (non-strange 
chiral condensate) decreases its value, causing the 
effective masses of the particles to decrease towards 
chiral symmetry restoration. The coupling constants for 
the baryons 36| are chosen to reproduce the vacuum 
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(5) 



with values of (/go- = g^^ = 4.0. At present we do not 
consider any possible quark- vectormeson couplings. 

A coupling of the quarks to the Polyakov loop is intro- 
duced in the thermal energy of the quarks. Their thermal 
contribution to the grand canonical potential fJ, can then 
be written as: 



and 



-T 



li 



^ (27r)3 



d^fcln 1 + (f>exp 



T 



(6) 



-T 



^ (27r)3 



d^fcln( l + ^^exp -^'*^^' ) (7) 



The sums run over all quark flavors, where 7^ is the 
corresponding degeneracy factor, E* the energy and [li 
the chemical potential of the quark. 
All thermodynamical quantities, energy density e, en- 
tropy density s as well as the densities of the different 
particle species p^, can be derived from the grand canon- 
ical potential. In our model it has the form: 



L 



V 



U 



(8) 



Here ^Ith includes the heat bath of hadronic and quark 
quasi particles. The effective potential [/(<&, $*, T) which 
controls the dynamics of the Polyakov-loop will be dis- 
cussed in the following. In our approach we adopt the 
ansatz proposed in [38|: 



U = -ia(r)$<I>* 



h{T)ln[l - 6$$* + 4($-^$*^) - 3($$ 



(9) 
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with a{T) = aoT* + fliToT^ + asT^T^, b{T) 



This choice of effective potential satisfies the Z(3) 
center symmetry of the pure gauge Lagrangian. In the 
confined phase, U has an absolute minimum at <i> = 0, 
while above the critical Temperature Tq (for pure gauge 
To = 270 MeV) its minimum is shifted to finite values 
of <&. The logarithmic term appears from the Haar 
measure of the group integration with respect to the 
SU(3) Polyakov loop matrix. The parameters ag, 01,02 
and 63 are fixed, as in [3§|, by demanding a first order 
phase transition in the pure gauge sector at Tq — 270 
MeV, and that the Stefan-Boltzmann limit is reached 
for T — > 00. Note that Tp remains a free parameter to 
adjust the actual critical temperature, of both phase 
transitions, when both, quarks and hadrons, couple to 
the scalar fields. 

As has been mentioned above, the Lagrangian of the 
chiral model contains dilaton terms to model the scale 
anomaly. These terms constrain the chiral condensate, 
if the dilaton is frozen at its ground state value xo- On 
the other hand, as deconfinement is realized, the expecta- 
tion value of the chiral condensate should vanish at some 
point. On account of this we will couple the Polyakov 
loop to the dilaton in the following way: 



X 



Xl (1 



(10) 



Hence, when the value of the Polyakov loop ap- 
proaches unity, the dilaton field will slowly vanish and 
allow the chiral condensate to also approach zero. 

Until now all hadrons are still present in the de- 
confined and chirally restored phase. Since we expect 
them to disappear, at least at some point above Tc, we 
have to include a mechanism that effectively suppresses 
the hadronic degrees of freedom, when deconfinement is 
achieved. 

In previous calculations baryons were suppressed by in- 
troducing a large baryon mass shift for non-vanishing $ 
0. 

In the following the suppression mechanism will be 
provided by excluded volume effects. It is well known 
that hadrons are no point-like particles, but have a finite 
volume. Including effects of finite-volume particles, in a 
thermodynamic model for hadronic matter, was proposed 
some time ago 41, 42, l43l 14^. We will use an ansatz 
similar to that used in 45. |46|. but modify it to also 



treat the point like quark degrees of freedom consistently. 

If one introduces a particle of radius r into a gas of the 
same particles, then the volume excluded is not just the 
simple spherical volume, but one-half times the volume 
of a sphere with radius 2r: 

v = \-^ir{2rf (11) 
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FIG. 1: (color online) The normalized order parameters for 
the chiral (red, solid line), and deconfinement (black, dashed 
line) phase transition as a function of T at ~ 0. The 
symbols denote lattice data from [s^l, using different lattice 
actions (asqdat and p4) and lattice spacings Nr. 



It is easy to understand that if all other particles also 
have a radius r then the excluded volume is much bigger 
than just the volume of a single particle. 
We expect the volume of a meson to be smaller than of 
a baryon we have to introduce the quantity Vi which is 
the volume excluded of a particle of species i. Since we 
only distinguish between hadronic baryons, mesons and 
quarks. Consequently Vi can only assume three values: 

VQuark = 
'^Baryon — ^ 
VMeson ^ v/a 

(12) 

where a is a number larger than one. In our calculations 
we assumed it to be a = 8, which would mean that the 
radius r of a meson is half that of a baryon. Note that 
we neglect any possible Lorentz contraction effects on the 
excluded volumes as introduced in [i^, [3 ■ 

The modified chemical potential Jli which is connected 
to the real chemical potential /i^, of the i-th particle 
species, is obtained by the following relation: 



Mi = Mi 



V, P 



(13) 



where P is the sum over all partial pressures. All thermo- 
dynamic quantities can then be calculated with respect 
to the temperature T and the new chemical potentials fli . 
To be thermodynamically consistent, all densities (e^, pi 
and s'i) have to be multiplied by a volume correction fac- 
tor /, which is the ratio of the total volume V and the 
reduced volume V' , not being occupied,: 



V 



(14) 
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FIG. 2: (color online) Total particle number densities for the 
different particle species as a function of T at /is = 0. The 
black line shows the total number of quarks+antiquarks per 
volume while the green (dotted) line refers to the total meson 
density and the red (dashed) line to the number density of 
hadronic baryons+antibaryons. 



The actual densities then are: 



= f P^ 

= E/ 



(15) 

(16) 
(17) 



Note that in this configuration the chemical potentials 
of the hadrons are decreased by the quarks, but not 
vice versa. In other words as the quarks start appearing 
they effectively suppress the hadrons by changing their 
chemical potential, while the quarks are only affected 
through the volume correction factor /. 
Our implementation of finite-volume corrections as out- 
lined above is a simple approach with as few parameters 
as possible and can be improved upon in various ways. 
For one, hadrons differ in size. The size of a hadron 
could even be density or temperature dependent [49j. In 
addition, the excluded-volume parameter of a particle 
does also depend on the density of the system (at dense 
packing a particle excludes effectively less volume). 
However, one should regard the variables v and a as 
effective parameters for capturing the qualitative effect 
of an excluded volume correction, which suppresses the 
hadrons in the quark phase. 



RESULTS AT VANISHING NET BARYON 
DENSITY 

In this paper we will concentrate on the properties of 
the model at vanishing chemical potential. Here different 



FIG. 3: The fraction of the total energy density that can be 
assigned to the quark-gluon phase {cqop contains the energy 
of the quarks and the Polyakov potential) as a function of T 
at /is = 0. 



lattice calculations suggest a crossover from the hadronic 
to the quark phase. Different lattice groups obtain differ- 
ent results for the phase transition t emp erature ranging 



from Tc = 160 MeV to 200 MeV [50|,l51|- For all follow- 
ing results we set Tq, the free parameter of the Polyakov- 
potential, to Tq = 235 MeV and the excluded volume 
parameter v = Ifm'^. This leads to a critical tempera- 
ture of Tc ~ 183 MeV {Tc is defined as the temperature 
with the largest change in the order parameters as a func- 
tion of the temperature) . 

The lattice data referred to in the following sections are 
taken from the HotQCD collaboration [H^. Here differ- 
ent actions (p4, asqtad) and lattice spacings {Nr = 6,8) 
were compared. Note that the transition region extracted 
from the lattice data lies between 185 and 195 MeV. 
Fig. [T] shows the temperature dependence for the or- 
der parameters of both, the deconfinement (<&), and chi- 
ral (a) phase transition, extracted from our model and 
compared to lattice data. Both order parameters change 
smoothly with temperature. The critical temperature is 
found to be equal for both phase transitions . Note that 
the value of the chiral condensate a only very slowly ap- 
proaches zero. This originates from the dilaton contribu- 
tion to the scalar potential in this model, which includes 
a repulsive term for small values of a. 

Fig. [3 shows the total densities of quarks plus anti- 
quarks (black solid line), mesons (green dotted line) and 
baryons plus antibaryons (dashed red line). Below the 
critical temperature hadrons are the dominant degree of 
freedom. When the quark number increases around Tc, 
they begin to suppress the hadrons. It is remarkable 
that the hadrons are still present, and not negligible, up 
to about 2.0 Tc [H^. Especially the mesons contribute 
strongly to all thermodynamic quantities, since they are 
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FIG. 4: (color online) Three times the pressure (red dashed 
line) and energy density (solid line) over as a function of 
T at hb = 0. The green dotted line indicates the Boltzmann 
limit for an ideal gas of three massless quarks and gluons. The 
symbols denote lattice data from fs^, using different lattice 
actions (asqdat and p4) and lattice spacings A'^^. 



FIG. 5: (color online) Energy density minus three times the 
Pressure over T* as a function of T at fj,B = 0. Also re- 
ferred to as the interaction measure. The dashed line depicts 
the interaction measure using the uncorrected energy density 
e. The symbols denote lattice data from ^53|, using different 
lattice actions (asqdat and p4) and lattice spacings A*',-. 



quite less suppressed than the baryons {vm < vb)- 
Above 2 Tc the hadrons are efFectively squeezed out of 
the system by the presence of the quarks. 

To emphasize this change in degrees of freedom, Fig. 
[3] shows the fraction of the total energy density which 
stems from the quarks and gluons (more precisely the 
Polyakov potential). As expected for a crossover both 
degrees of freedom (hadrons and quarks) are present in 
the temperature range from 0.75 — 2 Tc. Around Tc the 
fraction of the energy density, due to quarks and gluons 
increases rapidly. It converges to unity at around 2 times 
Tc. 

Let us now take a closer look at different thermo- 
dynamic quantities. Fig [?] displays the energy density 
(black curve) and three times the pressure (red dashed 
curve), both over T"* compared to lattice data [H^l- In 
the limit of infinite temperature, both quantities should 
converge to the Stefan-Boltzmann limit of an ideal gas 
of quarks and gluons. This limit is indicated as a green 
dashed line. The strong increase in energy density 
around Tc reflects the rapid change of the relevant 
degrees of freedom. At three times the critical tem- 
perature the energy density is slowly converging to the 
Stefan-Boltzmann limit, while the pressure is converging 
even slower as it was also observed in PNJL calculations 
54| . At temperatures below Tc our calculation gives 



larger values for the pressure and the energy density, 
which is not surprising as the correct description of the 
hadron contribution is problematic in the case of lattice 
calculations 



factor / of the excluded volume corrections. As has been 
shown above, the hadronic contribution to the densities 
disappears only at two time Tc and therefore they still 
exclude some portion of the volume for the quarks. The 
'dip' therefore indicates the disappearance of volume cor- 
rection factors for the quark phase. 

In the high temperature limit, where only the quarks 
(and gluons) remain in the system the energy density 
and pressure both slightly exceed the data from lattice 
calculations. 

Figure [H] displays the difference of the energy density 
and three times the pressure over T'^ (black solid line). 
This quantity is also referred to as the 'interaction mea- 
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Around 1.5 Tc one can observe a slight 'dip' in the en- 
ergy density. This 'dip' is connected to the correction 



FIG. 6: (color online) The speed of sound squared, as a func- 
tion of T/Tc at jUs = 0. 1/3 is the ideal gas limit. 
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FIG. 7: (color online) The Taylor coefficient C2 (red, solid FIG. 8: (color online) The Taylor coefficient C4 (red, solid 
line) as a function of T/Tc at ~ 0. The symbols denote line) as a function of T/Tc at /is = 0. The symbols denote 
lattice data from [S^). lattice data from [27} . 



sure' in lattice calculations. In the Stefan-Boltzmann 
limit it is 0, while it shows a peak slightly above Tc- 
The height of the peak in our model is comparable to 
the lower bound from lattice studies 5^, while its value 
at large T is a little bit above that from lattice calcula- 
tions, because chiral restoration is not fully achieved in 
our model. To point out the effect of the volume correc- 
tion factor on this quantity, Figure [5] also displays the 
uncorrected interaction measure (e — 3p) /T^ as a grey 
dashed line. 

An important property of a hot and dense nuclear 
medium is the speed of sound (cs): 



dp 

de 



(18) 



It is not only closely related to expansion dynamics but 
also controls the way perturbations (sound- and shock- 
waves) travel through the fireball [55] . Fig [6] shows the 
speed of sound squared as a function of temperature. As 
the temperature increases towards Tc one can clearly ob- 
serve a softening of the EoS due to the crossover. At very 
high temperature the speed of sound again converges to- 
ward its ideal gas limit of — > 1/3. The dip above is 
again related to the excluded volume corrections. Note 
that even though the change of degrees of freedom from 
hadrons to quarks proceeds as a crossover, there is still a 
substantial softening (i.e. goes down to 0.07)! 



EXTRAPOLATION TO FINITE BARYO 
CHEMICAL POTENTIALS 

Lattice results at finite chemical potentials (/is — 
3fiq 7^ 0) are often obtained as Taylor expansion of the 
thermodynamic quantities in the parameter fi/T around 



zero chemical potential and T = Tc- We perform the 
same expansion in the present model to compare with 
lattice results. In the Taylor expansion of the pressure 
p = —fl/V with respect to the coefficients c„ follow 
from: 



rp4 



CniT) 



E«..m(f)' 



n=0 



1 &\p{T,^Jiq)/T^) 



n! ^{^iJTy 



(19) 
(20) 



Because the derivatives of the pressure with respect to 
/ig/T are difficult to be obtained directly from a mean 
field approach (see [11]), we calculate the pressure as 
function of T and /ig explicitly and then extract the 
derivatives around ^iq = 0. 



The results for the second moment are shown in Fig 
[7] and compared to lattice calculations taken from 27 ) . 
The second order coefficient C2 , exhibits an enhancement 
below Tc, due to the hadronic contributions and a 
minimum just above Tc- This is not observed in lattice 
calculations, because the origin of this dip is the presence 
of the additional hadronic degrees of freedom and the 
excluded volume effects in the hadronic part of the 
model. Thus it is not surprising to see a discrepancy to 
the lattice data. 

The fourth order coefficient C4 is shown in Fig [H As 
expected, a peak is observed around Tc- Again the 
absolute height of the peak exceeds that from lattice 
calculations. This, again, stems from the fact that 
our model includes many hadronic degrees of freedom 
around Tc- The minimum above Tc is also visible in the 
forth order coefficient. 

Note that the quark number susceptibility is also closely 
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related to the Taylor coefficients and therefore shows a 
clear maximum at Tc- 

As our model does not sustain the difficulties of lattice 
calculations, when going to finite net baryon densities, it 
is straight forward to extend our investigations to finite 
chemical potentials {^b = 7^ 0). These calculations 
at finite densities are in progress and will be subject of 
future publications. 



CONCLUSION 

We present a novel approach for modeling an EoS that 
respects the symmetries underlying QCD, and includes 
the correct asymptotic degrees of freedom, i.e. quarks 
and gluons at high temperature and hadrons in the low- 
temperature limit. We achieve this by including quarks 
degrees of freedom and the thermal contribution of the 
Polyakov loop in a hadronic chiral sigma-omega model. 
The hadrons are suppressed at high densities by excluded 
volume corrections. Nevertheless, we observe a substan- 
tial hadronic contribution to the EoS up to temperatures 
of 2 times Tc. 

We can show that the properties of the EoS are in qual- 
itative agreement with lattice data at fis — 0. Various 
quantities, like the pressure and energy density, are in 
good agreement with lattice data. Deviations from lat- 
tice results , for example the Taylor coefficients of the 
expansion in fJ./T, can be explained by the hadronic con- 
tributions and volume corrections. In spite of a contin- 
uous phase transition, we obtain a considerably smaller 
value for the speed of sound around Tc (c^ ~ 0.072) when 
compared to lattice calculations 
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